Non-perturbative calculation of the two-loop Lamb shift in Li-like ions 
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A calculation valid to all orders in the nuclear- strength parameter is presented for the two-loop Lamb shift, 
notably for the two-loop self-energy correction, to the 2p -2s transition energies in heavy Li-like ions. The 
calculation removes the largest theoretical uncertainty for these transitions and yields the first experimental 
identification of two-loop QED effects in the region of the strong binding field. 
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The Coulomb field of heavy nuclei provides a unique op- 
portunity for testing the strong-field regime of bound-state 
quantum electrodynamics (QED). The most obvious candi- 
date for such a test is an H-like ion, whose theoretical descrip- 
tion is the simplest one. Measurements of the Is Lamb shift 
in H-like uranium, the heaviest naturally occurring element, 
have recently achieved an accuracy of 4.6 eV [1]. This cor- 
responds to a fractional accuracy of L7% with respect to the 
total Is QED contribution. Such measurements yield a test of 
bound- state QED at the one-loop level (i.e., to first order in 
the fine structure constant a), but they are yet insensitive to 
the two-loop QED corrections, which are of primary theoreti- 
cal interest at present. 

By contrast, measurements of the 2p-2s transition energies 
in heavy Li-like ions l^l, Sll, B] have lately reached a frac- 
tional accuracy of 0.03% with respect to the total QED con- 
tribution. This corresponds to a 10% sensitivity of the ex- 
perimental results to the two-loop QED contribution. These 
measurements provide an excellent possibility for identifica- 
tion of the two-loop Lamb shift and for testing the bound- state 
QED up to second order in a in the strong-field regime. 

The theoretical description of Li-like ions is complicated 
by the presence of the electron-electron interaction. For heavy 
ions, this interaction can be successfully accounted for within 
the perturbative expansion in a small parameter 1/Z (Z is the 
nuclear charge number). By calculating a few terms of this ex- 
pansion, one can rigorously describe the electron correlation 
and the screening of one-loop QED corrections with the ac- 
curacy sufficient for identification of the two-loop Lamb shift. 
Such a project has recently been accomplished in [ 6, 7] . Based 
on these calculations, an "experimental" value of —0.23 eV 
was inferred in |5] for the 2s two-loop Lamb shift in Li-like 
uranium. A similar determination of the two-loop Lamb shift 
was earlier presented in |7] (based on the measurement 0]) 
for the 2^3/2-25 transition energy in bismuth. 

The subject of our present interest is the set of two-loop 
one-electron QED corrections (also referred to as the two-loop 
Lamb shift), graphically represented in Fig. [T] These cor- 
rections have lately been extensively investigated within the 
perturbative expansion in the nuclear- strength parameter Za 
S S [13, [m [12] • Such studies, however, do not provide re- 
liable information about the magnitude of the two-loop Lamb 
shift in heavy ions like uranium, where the parameter Za ap- 
proaches unity. Our present investigation is addressed primar- 



ily to high-Z ions and will be performed non-perturbatively, 
i.e., without an expansion in Za. The only exception will be 
made for the diagrams in Fig.[nh)-(k), for which we will ex- 
pand the fermion loops in terms of the binding potential. We 
will keep the leading term of the expansion and refer to this 
as the free-loop approximation. In the one-loop case, such ap- 
proximation corresponds to the Uehling potential and yields 
the dominant contribution even for high-Z ions like uranium. 

Necessity for a non-perturbative calculation of the two- 
loop Lamb shift became clear already in the beginning of 
the 1990s, after the famous measurement 10] of the 2pi/2-25 
transition energy in U^^+ with an accuracy of 0.1 eV. Quite 
soon, calculations of all diagrams in Fig.[nd)-(k) were accom- 
plished [il3|] [although the graphs (i) and (k) were calculated 
in the free-loop approximation only]. Three remaining dia- 
grams (a)-(c), referred to as the two-loop self-energy correc- 
tion, turned out to be much more difficult to evaluate. The cal- 
culation for the Is state extended over a decade H Hill], 
with the first complete evaluation presented in uh and later 
extended in [18]. 

In this Letter, we present our calculation of the two-loop 
self-energy correction for the 2s, 2pi/2, and 2^3/2 states of 
several high-Z ions. We also evaluate the remaining diagrams 
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FIG. 1: Two-loop one-electron QED corrections. Gauge-invariant 
subsets are referred to as SESE (a-c), SEVP (d-f), VPVP (g-i), 
S(VP)E (k). 
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TABLE I: The two-loop self-energy correction, in terms of F{Za). 



z 


Is 


2s 


2pi/2 


2p3/2 


60 


-1.666(19) 


-1.976(70) 


0.222(72) 


0.02 (14) 


70 


-1.923(18) 


-2.453 (60) 


0.192 (60) 


-0.082 (93) 


83 


-2.360(15) 


-3.296 (38) 


0.133 (38) 


-0.175 (66) 


92 


-2.806(12) 


-4.218(34) 


0.012(32) 


-0.241 (48) 


100 


-3.392(14) 


-5.455 (68) 


-0.214(32) 


-0.282 (52) 



in Fig. [Hand obtain results for the total two-loop Lamb shift, 
this being previously the largest uncalculated correction for 
the 2p-2s transition energies in heavy Li-like ions. Our cal- 
culation significantly refines theoretical predictions for these 
transition energies and provides a test of bound- state QED 
theory in the strong-field regime up to second order in a. 

Another important aspect of our calculations is associated 
with their implications for the hydrogen Lamb shift. In the 
previous evaluation of the Is two-loop self-energy correction 
LI 8.1 . we demonstrated that an extrapolation of our results to 
Z = 1 yields a value that disagrees with the analytical result 
to order a^(Za)^ [11]. In view of this disagreement, it will be 
of interest to compare our non-perturbative values with the an- 
alytical result to order {Zotf for the normalized difference 
of the 2s and \s energy shifts, = 8(^£^2s — ^^Xs- Such 
a comparison can provide us with new information about the 
discrepancy for the \s state, since, within the perturbative ap- 
proach, the difference is understood much better than the 
energy shift of a single ns state [19]. 

The calculational scheme for the evaluation of the two-loop 
self-energy correction was developed for the ground state in 
our previous studies Ij/Zi [l8|]. With this work, we extend 
it to the excited states. The general procedure for isolation 
and cancelation of divergences is similar to that for the \s 
state, but the actual calculational scheme requires substantial 
modifications due to a more complicated pole and angular- 
momentum structure of expressions involved. Details of our 
calculation will be published elsewhere; in this Letter we con- 
centrate on presentation and analysis of the results obtained. 

The two-loop self-energy correction to the energy is conve- 
niently represented in terms of the function F defined by 

5S = mc2n' (^F(Za), (1) 

where n is the principal quantum number. The numerical re- 
sults obtained for the n = 1 and 2 states are listed in Ta- 
ble H The calculation was performed for the point model of 
the nuclear-charge distribution. 

As an intermediate step in our calculation, we had to con- 
sider the irreducible part of the diagram in Fig.IHa), also de- 
noted as the loop-after-loop correction, previously calculated 
in fl3]. We report a good agreement with the previous re- 
sults for the \s state but find a discrepancy for 2s and 2pi/2 
states. For Z = 92, we obtain -0.090 eV and -0.030 eV 
for the 2s and the 2pi/2 state, respectively, which should be 
compared with the values of —0.069 eV and 0.014 eV from 
|[l4ll . respectively. The apparent reason for this disagreement 
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FIG. 2: Non-perturbative results for the higher-order remainder 
G^'^'{Z) for the differences As (squares) and Ap (dots) and their 
limiting values at Z = denoted by the cross for As and by the 
diamond for Ap. 

is a sign error in the contribution due to the imaginary part of 
the self-energy operator. 

The consistency of our numerical values can be tested by 
comparing them with analytical results obtained within the 
perturbative approach. The Za expansion of the function F 
reads 

F{Za) = B40 + {Za) B^o + {Zaf [L^ Bqs 

+ Be2 ^LBei^ G^-"" (Za)] , (2) 

where L = ln[(Za)~^] and G^°- is the remainder, 
G^°(Za) = Beo^ Za{---). The coefficients B^q-Bqi are 
presently known for all states of our interest |8lfi[i3] • The co- 
efficient was calculated for the specific differences of en- 
ergy shifts As = S 8E2S - SEis and A^ = ^^2^3/2 - ^^2^1/2 
ifTol [T2h . A calculation of the dominant part of Bqo{1s) and 
(25) was reported in |11], together with an estimate of 
unevaluated contributions to this order. 

We would like now to isolate the contribution of the higher- 
order remainder G^°" from our numerical results. Obviously, 
such isolation leads to a significant loss of precision, which 
grows fast when Z decreases. Moreover, the data presented in 
Table U for the function F are already a result of a significant 
(and Z-dependent) cancelation, since individual contributions 
to the energy do not exhibit the physical Z^ dependence but 
scale typically as Z or Z^ . This indicates that an analysis of 
the Z behavior of the higher-order remainder G^ °" provides a 
stringent test of the consistency of the numerical results. 

In Fig. [2I we present our results for the higher-order re- 
mainder G^ ° (Z) for the differences A5 and Ap, together 
with their limiting values G^ °^) = Bqq, which are 14.1(4) 
for As and —0.3611 for Ap |10', The comparison pre- 
sented in the figure indicates that our non-perturbative results 
are in good qualitative agreement with the perturbative expan- 
sion coefficients known for the differences As and Ap. A 
quantitative comparison is presently impossible due to a lack 
of numerical data in the low-Z region. To perform a non- 
perturbative calculation for excited states of low-Z ions is a 
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TABLE II: Individual two-loop contributions to transition energies 
in Li-like bismuth and uranium, in eV. 



Subset 


Fig. 


2p3/2-2s, Z = 83 


2pi/2-2s, Z = 92 


SESE 


(a)-(c) 


0.145 (4) 


0.296 (3) 


SEVP 


(d)-(f) 


-0.095 


-0.187 


VP VP 


(g) 


0.016 


0.035 


VP VP 


(h),(i) 


0.067(25) 


0.101(46) 


S(VP)E 


(k) 


-0.012(24) 


-0.022 (45) 


Total 




0.120 (35) 


0.223 (64) 



difficult problem, whose solution apparently requires devel- 
opment of new calculational technique. For the Is state, how- 
ever, such a calculation is less problematic and we were able to 
carry out a direct evaluation for Z as low as 10 in our previous 
investigation [18]. Extrapolating our Is results to Z = 1 in 
that work, we found a disagreement with the analytical value 
for the coefficient Bqq{1s) |11]. Our present calculation for 
the difference As and the agreement observed with the ana- 
lytical expansion coefficients in this case can be considered as 
an evidence in favor of reliability of our previous results for 
the Is state. 

In Table nil we present the results of our calculations of all 
two-loop corrections depicted in Fig. [T] for the 2^3/2-25 tran- 
sition in Bi^^+ and the 2pi/2-25 transition in U^^+, for which 
most accurate experimental data are available. The contribu- 
tion of the SESE subset is taken from Table Jl The diagrams 
in Fig.[nd)-(g) were calculated rigorously to all orders in Za, 
whereas the diagrams in Fig. [Ilh)-(k) were evaluated within 
the free-loop approximation, i.e., keeping the first nonvan- 
ishing contribution in the expansion of the fermion loops in 
terms of the binding potential. The error bars specified for 
these corrections are estimations of uncertainty due to the ap- 
proximation employed. They were obtained by multiplying 
the contribution of Fig. \V[h,i) by a factor of (Za)'^ and that 
of Fig. mk), by a factor of 3 (Za). The factor of 3 {Za) 
in the latter estimation arises as a ratio of the leading-order 
contribution beyond the free-loop approximation for the di- 
agram (k), -0.386 (a/7r)2(Za)^ ||20|], and the leading-order 
contribution within this approximation, 0.142 (a/7r)^(Za)^ 
1I2TI1 . The finite nuclear size effect was taken into account 
in our evaluation of the diagrams in Fig.[nd)-(i), whereas the 
other diagrams were calculated for the point nuclear model. In 
the case of uranium, our results for the diagrams with closed 
fermion loops are in good agreement with those reported pre- 
viously 1 22]. 

We now turn to the experimental consequences of our cal- 
culations. In Table [nil we collect all available theoretical con- 
tributions to the the 2^3/2 -25 transition energy in Bi^^+ and 
to the 2pi/2-25 transition energy in U^^+. The entry labeled 
"Dirac value" represents the transition energies as obtained 
from the Dirac equation with the Fermi-like nuclear potential 
and the nuclear-charge root-mean- square (rms) radius fixed as 
<r^^/2= 5.851(7) Em for uranium and 5.521(3) Em for bis- 
muth 1I23I] . The dependence of the Dirac value on the nuclear 
model was conservatively estimated by comparing the results 
obtained within the Fermi and the homogeneously-charged- 



TABLE III: Various contributions to transition energies in Li-like bis- 
muth and uranium, in eV. 





2p3/2-2s, Z = 83 


2pi/2-2s, Z = 92 


Dirac value 


2792.21 (3) 


-33.27(9) 


One-photon exchange 


23.82 


368.83 


Two-photon exchange 


-1.61 


-13.37 


Three-photon exchange 


-0.02(2) 


0.15(7) 


One-loop QED 


-27.48 


-42.93 


Screened QED 


1.15(4) 


1.16(3) 


Two-loop QED 


0.12(4) 


0.22 (6) 


Recoil 


-0.07 


-0.07 


Nuclear polarization 




0.04(2) 


Total theory 


2788.12 (7) 


280.76 (14) 


Experiment 


2788.14(4) [3] 


280.645 (15) 1^ 



sphere models (2^. We have checked that a wide class of 
more general models for the nuclear-charge distribution yields 
results well within the error bars obtained in this way. 

The next 3 lines contain the corrections due to the one-, 
two-, and three-photon exchange, respectively. QED values 
for the two-photon exchange correction were taken from our 
previous evaluations |6, 25]. The results for the three-photon 
exchange correction were obtained in this work within many- 
body perturbation theory (MBPT), with retaining the Breit in- 
teraction to the first order only. For uranium, we report good 
agreement with the previous evaluations of this effect ll26n . 
The error ascribed to the three-photon exchange correction is 
due to incompleteness of the MBPT treatment. It was esti- 
mated by calculating the third-order MBPT contribution with 
two and more Breit interactions for each state involved in the 
transition, adding these contributions quadratic ally, and mul- 
tiplying the result by a conservative factor of 2. 

The entry labeled "One-loop QED" represents the sum of 
the first-order self-energy and vacuum-polarization correc- 
tions calculated on hydrogenic wave functions ll22h . The next 
line ("Screened QED") contains the results for the screened 
self-energy and vacuum-polarization corrections ||27|] . The 
uncertainty ascribed to this entry is the estimation of higher- 
order screening effects; it was obtained by multiplying the 
correction by the ratio of the entries "Screened QED" and 
"One-loop QED". The entry "Two-loop QED" contains the 
results for the two-loop Lamb shift obtained in the present 
investigation. The next two lines contain the values for the 
relativistic recoil correction ll28h and the nuclear polarization 
correction [29]. 

Table |llll shows that now, after our calculation of the two- 
loop Lamb shift, the total theoretical uncertainty is signifi- 
cantly influenced by the error of the finite-nuclear- size effect. 
It should be noted that a certain concern exists in the commu- 
nity about the accuracy of the theoretical description of this 
correction. In particular, it was pointed out 1 30] that little is 
known about systematical effects in experimental determina- 
tion of nuclear rms radii. In the absence of detailed investiga- 
tions of such effects, we consider the errors of rms radii ob- 
tained in [23] by averaging all experimental results available 
(including both muonic-ions and electron- scattering data) to 
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be presently the most reliable estimates and employ these val- 
ues in our calculations. 

The comparison presented in Table [Till demonstrates that 
our total results agree well within the error bars specified with 
the experimental data for bismuth and uranium. The theoreti- 
cal accuracy is significantly better in the former case, which is 
the consequence of the fact that the finite nuclear size effect is 
smaller and the nuclear radius is known better. Our result for 
the 2p3/2-2s transition in bismuth can also be compared with 
the value of 2787.96 eV obtained by Sapirstein and Cheng | 7]. 
The difference of 0. 16 eV between the results is mainly due to 
the two-loop Lamb shift contribution (0.12 eV) which is not 
accounted for in [7] . 

We conclude that inclusion of the two-loop Lamb shift is 
necessary for adequate interpretation of the experimental re- 



sult in the case of bismuth, whereas for uranium the two-loop 
Lamb shift is significantly screened by the uncertainty due to 
the nuclear charge distribution. Comparison of the theoretical 
and experimental results for bismuth yields the first identifica- 
tion of the two-loop QED effects in the region of strong bind- 
ing field, which is the first step toward the test of the strong- 
field regime of bound- state QED at the two-loop level. 
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